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Abstract: 

Coupled discrete models abound in several areas of physics. Here we provide an extensive set of 
exact quasiperiodic solutions of a number of coupled discrete models in terms of Lame polynomials of 
arbitrary order. The models discussed are (i) coupled Salerno model, (ii) coupled Ablowitz-Ladik model, 
(iii) coupled cp'^ model, and (iv) coupled (f)^ model. In all these cases we show that the coefficients of 
the Lame polynomials are such that the Lame polynomials can be reexpressed in terms of Chebyshev 
polynomials of the relevant Jacobi elliptic function. 
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1 Introduction 



In a recent paper we have obtained solutions of a number of coupled discrete models in terms of Lame 
polynomials of order one and two. The purpose of the present paper is to show that in the same models 
one can in fact obtain solutions in terms of Lame polynomials (and hence trigonometric and hyperbolic 
polynomials) of arbitrary order. In particular, we obtain solutions of (i) Coupled Salerno model (ii) coupled 
Ablowitz-Ladik (AL) model (iii) coupled (j)^ model (iv) coupled (f)^ model, in terms of Lame (and hence 
trigonometric and hyperbolic) polynomials of arbitrary order. As an illustration, we confine our discussion 
to coupled Salerno model and then show how similar solutions also exist in coupled Ablowitz-Ladik, coupled 
(f)^ and coupled (j)^ models. Quite remarkably, we find that the coefficients of the Lame polynomials are such 
that the Lame polynomials can be reexpressed as Chebyshev polynomials of the relevant Jacobi elliptic 
function. 

There are many physical situations where a discrete field theory is appropriate to model the phenomena 
of interest with a specific coupling between the two fields. This fact serves as one of the main motivation 
for the results presented here. A technologically important phenomenon of current intense interest is 
the coexistence of magnetism and ferroelectricity (i.e., magnetoelectricity) in a given material. This is a 
highly desired functionality in applications involving cross-field response such as electric field control of 
magnetism and magnetic field control of electric polarization, switching and actuation. More generally, 
presence of two or more such ferroic properties (e.g. polarization, magnetization, strain) is referred to as 
multiferroic behavior j^. In recent years, two different classes of (single phase) multiferroics, namely the 
orthorhombically distorted perovskites [5] and rare earth hexagonal structures [6], have emerged in addition 
to other related structures. The hexagonal materials show axial polarization and magnetic domain walls 
in the basal planes. The latter can be modeled by a coupled (f)'^ model [7] in the presence of a magnetic 
field. Coupled cj)'^ models [H [9l [10] are useful in the study of many ferroelectric and other second order 
phase transitions. For multiferroic materials the relevant coupled (j)'^ model [7j has a biquadratic coupling. 
In contrast, the coupled (f)^ model for a surface phase transition with hydration forces jlOj . which arises 
in the biophysics context, has a bilinear coupling. Other types of couplings are also known for structural 
phase transitions involving strain tensor components |llj . 
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There are many examples of coupled discrete Ablowitz-Ladik, coupled discrete Salerno and coupled 
saturated discrete nonlinear Schrodinger (DNLS) models known in the literature [12 \ \1'6 \ [H]. In addition, 
there are analogous coupled models known in the field theoretic contexts [15[ [TC] . Several related models 
have been discussed in the literature and their soliton solutions have been obtained [TTJ [TSJ [TOl EH ED ESI 
E31 ED including periodic ones ESI EZ]- 

The paper is organized as follows. In Section II we provide the solutions for the coupled Salerno model 
in terms of Lame polynomials of order three and four. Based on these results as well those obtained in pp, in 
Section III we generalize these results and conjecture solutions in terms of Lame (and hence trigonometric 
and hyperbolic) polynomials of arbitrary order. In Section IV we show that these Lame polynomials can 
be reexpressed as Chebyshev polynomials of the relevant Jacobi elliptic function. This also proves that our 
proposed Lame polynomials of arbitrary order are indeed solutions of the coupled equations. In Section V 
we show how the coupled Ablowitz-Ladik, coupled discrete (j)^, and coupled (/>^ models also admit solutions 
in terms of Lame polynomials of arbitrary order. Section VI contains the summary of main results and 
possible future directions. 



2 Coupled Salerno Model 



As discussed in a recent paper [T], the field equations of the coupled Salerno model are given by 



idUn/dt + [Un+l + - 2Un] + (/Ul I Un |^ +1^2 \ Vn P) 



vi — 2^1 

Un+l + Un-l H 



0, (1) 



idvn/dt + 



Vn+l + Vn^l - 2 + 



+ I Un P +/i2 I Vn ?) 



V2 — 2/^2 

■Wn+l + Vn-\ H Vr. 

ii2 



= 0. 

(2) 



In (HI) we have already obtained solutions of these coupled equations in terms of Lame polynomials of 
order one and two. We now show that the same model also admits Lame polynomial solutions of arbitrary 
order. 

As in [1], we start with the ansatz 



Un = fn exp [-i{uJit + 6i)] , Vn= Qn exp [-i{u2t + 82)] , 



(3) 
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then it is easily shown that the above coupled equations take the form 

(Wi - 2)/„ + + /„_l)(w/n + + 1) + (m/n + f^2gl)^^^^^^ fn = , (4) 
U2-2 + - — ) Qn + {^llfl + fl2gl + l)(fi'n+l + Qn-l) + {^llfn + f^2gn)— = . (5) 

A* 1 y 

It is clear from Eqs. ^ and ([5]) that in general these coupled equations will have exact solutions if 

l + /^l/n+^25n = 0) (6) 

and in that case 

OJl = , UJ2 = 2~ • ( 

We now show that Eq. ^ has solutions in terms of Lame polynomials of arbitrary order. In particular, 
we show that at every order there are three distinct solutions satisfying Eq. As a first step, let us 

explicitly obtain solutions of Eq. ([6]) in terms of Lame polynomials of order three and four. (It may be 
noted that we have already obtained such solutions before in terms of Lame polynomials of order one and 
two d].) 

2.1 Lame Polynomial Solutions of Order Three 

It is easily checked that one of the exact solution to Eq. ([6]) is 

fn = dn[/3(n + C2),m](^dn2[/3(n + C2),m] + 5) , Qn = \/^sn[/3(n + C2), m](Cdn2[/3(n + C2), m] + L>) , (8) 
provided 

/Zi,/X2<0, \^li\A^ = \^l2\C'^ , |/X2|Z)^ = 1, 

+ = 1, \^xi\{A^ + 2AB) = 2CD\^l2\. (9) 

Here C2 is an arbitrary constant signifying discrete translation invariance. Also notice that for this solution, 
the width (3 is also completely arbitrary. We shall see that all the solutions discussed in this paper are 
valid for arbitrary C2 and width /3. On solving, we find that 

^\A = A,B = --^A, ^\C = A, C = -AD. (10) 



Since the field Eqs. (j3|) and ([5|) are invariant under — >• ±/„, gn — )• ±5n) one can trivially write three other 
solutions from here. Note that /3 is completely arbitrary. Using the fact that dn(x, m) has period 2K{m) 
while cn(x,m) and sn(x,m) are periodic functions with period 4:K{m), it follows that for the solution ([5]), 
fn,9n and hence Un,Vn satisfy the boundary condition 

U 2K(m) = Un , f 4K(m) = ■ (H) 

Here K{m) is the complete elliptic integral of the first kind. 
Another solution to Eq. ([6]) is 

fn = cn[/3(n + C2),m]{Acn^[P{n + C2),m] + B) , = sn[/3(n + cg), m](Ccn^[/3(n + C2),m] + D) , (12) 

provided A,B,C,D are again given by Eq. (llOp . Note however that the solution (I12p is distinct from the 
solution dS]). In particular, while for the solution ([5]), Un,Vn satisfy the boundary condition (jlip . for the 
solution fn,gn and hence Un,Vn satisfy the boundary condition 

, 2K(m) = Un , t'^ , 2K(m) = fn • (13) 

However, in the limit m = 1, both the solutions dH) and (|12p go over to the hyperbolic solution 

/„ = Asech^lPin + cs)] + Bsech[/3(n + C2)] , 5n = tanh[/3(n + C2)](Csech2[/3(n + C2)] + D) . (14) 

Further, in the limit m = 0, the solution ()12p goes over to the trigonometric solution 

/n = cos[/3(n + C2)](Acos2[/3(n + C2)] + 5) , gn = sin[/3(n + C2)](C cos2[/3(n + C2)] + £>) . (15) 

All the four solutions discussed above are valid only if /Ui,/i2 < 0. If instead ^i,/i2 have opposite signs 
then also one has solutions to Eq. ([6]). One such solution is 

J, 1 / A ^ \ ^ _ V^sn[/3(n + C2),m] / C 

dn[/3(n + C2),m] Vdn^[/3(n + C2), m] / ' dn[/3(n + C2), m] Vdn^[/3(^^ + C2), m] 

(16) 

provided < 0, /^2 > while A, B, C, D are still given by Eq. (jlOp . Note that for this solution /„, (7„ and 
hence f„ satisfy the boundary condition (jlip . Further, if we interchange and gn, then > 0, /.f2 < 0. 
In the limit m = 1, this solution goes over to the hyperbolic solution 

/„ = Acosh^lPin + C2)] + Bcosh[P{n + C2)] , 5n = sinh[/3(n + C2)](C cosh2[/3(n + C2)] + -D) . (17) 
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2.2 Lame Polynomial Solutions of Order Four 

It is easily checked that one of the exact sokition to Eq. ([6]) is 

/„ = Adn^[f3{n + cs), m] + Bdn^[P{n + ca), m] + C , 

On = \/^sn[/3(n + C2),m]dn[/3(n + C2),m](Z)dn2[/3(n + C2),m] +£;), (18) 

provided 

^ii,^2<0, l^nl^^ = |/U2|^^ \^Li\C^ = l, A{A + 2B)\^i\ = 2DE\^i2\, 

\^Xl\{A + Bf + 2AC\iii\ = \ii2\E'^ , 2BC\iii\ = -\^12\E'^ . (19) 

On solving, we find that 

yi^C = 1,A = = 8C, vlAd^ = -4, D = -2E. (20) 

For the solution p^ . /n,5n and hence Un,Vn satisfy the boundary condition pT]) . 
Another solution to Eq. ([6]) is 

/„ = Acn'^[/3{n + C2),m] + Bcn^[l3{n + C2),m\ + C , 

Qn = sn[/3(n + C2),m]cn[/3(n + C2), m](Dcn^[/3(n + C2),m] + iiJ) , (21) 

provided A,B,C,D,E are again given by Eq. (|2Up . Note however that the solution (|21|) is distinct from 
the solution (jlSp . In particular, while for the solution (llSp . Un,Vn satisfy the boundary condition (llip . for 
the solution (I2ip . fn,gn and hence Un,Vn satisfy the boundary condition (I13p . 

However, in the limit m = 1, both the solutions (jl8p and (j2ip go over to the hyperbolic solution 

= Asech^[/3{n + C2)] + 5sech2[/3(n + ca)] + C , 
gn = tanh[/3(n + C2)]sech[/3(n + C2)](Dsech2[/3(n + cs)] + E) , (22) 

Further, in the limit m = 0, the solution (|2ip goes over to the trigonometric solution 

fn = A cos^[/3(n + C2)] + S cos2[/3(n + C2)] + C , 

5„ = sin[/3(n + ca)] cos[/3(n + C2)](D cos2[/3(n + cs)] + E) , (23) 
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All the four solutions discussed above are valid only if fii, fi2 < 0. If instead Hi, IJ-2 have opposite signs 
then also one has solutions to Eq. One such solution is 



dn^[(3{n + C2),m] dn^[/3(n + C2), m] 



gn = V^sn[(3{n + C2),m]( ^ - + ^ ) , (24) 

Vdn [p[n + C2),m\ dn [p[n + C2),m\J 

provided < 0, > while A, B, C, D, E are still given by Eq. ()20p . For this solution, fn,gn and hence 
Un,Vn satisfy the boundary condition (llip . Note that if we interchange /„ and gn, then > 0,/i2 < 0. In 
the limit m = 1, this solution goes over to the hyperbolic solution 

/„ = Acosh^[/3{n + C2)] + 5cosh2[/3(n + C2)] + C , 

5„ = sinh[/3(n + C2)](£'cosh3[/3(n + C2)] +^cosh[/3(n + C2)]). (25) 



3 General Results 

Looking at the structure of the solutions in terms of Lame polynomials of order one to four, it is easy to 
generalize and write down the solutions of Eq. ([6]) in terms of Lame polynomials of arbitrary order. For 
this we need to divide the discussion into two parts depending on if we are considering Lame polynomials 
of odd or even order. 

Case I: Lame Polynomials of odd order 

One of the solutions can be written in the form (note n is an odd integer) 

(n+l)/2 

fn= Yl Akidn[/3{n + C2),m]f''-\ 

k=l 

(n+l)/2 

5n = V^sn[/3(n + C2),m] Bkidn[f3{n + C2),m]f^-\ (26) 

k=l 

For this solution, /„, gn and hence Un,Vn satisfy the boundary condition Note that there are (n+ 1) /2 

number of terms in both /„ and gn- 
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Another solution is given by 

(n+l)/2 

fn= Yl Mcr^[(3{n + C2),m]) 



2k— 1 
k=l 



(n+l)/2 

gn = sn[/3{n + C2),m] ^ 5fc(cn[/3(n + ca), m])^*^-^ . (27) 

k=l 

For this solution, fn,gn and hence Un,Vn satisfy the boundary condition 

In the limit m = 1, both the solutions ()26p and ()27p go over to the hyperbolic solution 

(n+l)/2 

/„= J] ^fc(sech[/3(n + C2)])2'=-i, 

k=l 

(n+l)/2 

ff„ = tanh[/3(n + C2)] Sfc(sech[/5(n + cs)])'^-^ . (28) 

k=l 

Further, in the limit m = 0, the solution ()27p goes over to the trigonometric solution 

(ri+l)/2 

/n= 5^ ^fc(cos[/3(n + C2)])2^-i 
fc=l 



(n+l)/2 

5n = sin[/3(n + C2)] 5fc(cos[/3(n + cs)])''^-^ (29) 

k=l 

Note that these four solutions are valid in case /^i,/i2 < 0. However, if fii and /i2 have opposite signs, 
say /ii < 0, /i2 > 0, then the solution is given by 

(n+l)/2 



2A:-1 



^ (dn[/3(n + C2),m]) 

(n+i)/2 ^ 

ffn = V^sn[l3{n + C2), m] V — — r^g^ . (30) 

^ (dn[p(n + C2), mj)^'^ ^ 



For this solution, fn,gn and hence Un,Vn satisfy the boundary condition (|lip . In the limit 7n = 1, this 
solution goes over to the hyperbolic solution 

(n+l)/2 

fn= Yl M^Oshmn + C2)]f^-\ 
k=l 

{n+l)/2 

(7„ = tanh[/3(n + C2)] Bfc(cosh[;g(n + 02)])'^"^ (31) 

fe=i 

On substituting any of the expressions for fn,gn as given by Eqs. (I26p to (I3ip in Eq. ([6]), we obtain 



n + 1 equations which determine the n + 1 parameters Ak^Bk- While all A^s are numbers in units of i — 

VImiI ■ 



all BkS are numbers in units of / For simplicity from now onwards we will merely give the numerical 
values of A^, and it is understood that they are in units of / and / , respectively. Some of the 
relations are 

^(n+l)/2 - ^{n+l)/2 ' 

2-4(„+i)/2^(„-l)/2 + Bf^+l)/2 = 2-B(n+l)/2-B(n-l)/2 , 

^(n-l)/2 + 2^{n+l)/2^(r!,-3)/2 + 2S(„+l)/2^(r!,-l)/2 = -^(n-l)/2 + 2S(r!,+l)/2-B(n-3)/2 > 
-^{n-l)/2 + 2^{n+l)/2^{n-5)/2 + 2-B(„+i)/2S(„_3)/2 + 2^(„_i)/2^(„-3)/2 
= 2-B(^_i)/2i?(n-3)/2 + 2-B{n+l)/2-B(,„_5)/2 , 

Bl = l, A\ + 2BiB2 = Bl, 2AiA2 + B^ + 2BiB3 = 2BiB2, 

Al + 2A1A3 - bI + 2B2B3 + 251^4 - 251^3 = . (32) 

On comparing these results with the exact expressions for n = 1,3, we conjecture the following general 
results for arbitrary odd n: 

^(n+l)/2 = -B(„+i)/2 = 2"~\ ^(n-l)/2 = -"-2""^ , -B(„_i)/2 = - - 2)2"~^ , 

_ n(n^-3) 5 (ra-3)(n-4) 5 

^{n-3)/2 - ^ , -D(n-3)/2 - 

n(n-4)(n-5) y (n - 4)(n - 5)(n - 6) 7 

^(n-5)/2 = ^1 2 , -0(„_5)/2 = 2 , 

(n+l)/2 (n+l)/2 

^ ^i = l, ^ Bi = n. (33) 

j=l i=l 

Further, depending on if n = 4A; + 3 or 4A; + 1 we have the following results: 

Ai = -n, Bi = -1, if n = ik + 3, A; = 0, 1,2, ... (34) 

Ai = +n, Bi = +1, if n = 4:k + l, A; = 0, 1, 2, ... . (35) 

In fact looking at these general results, we conjecture the following expressions for various Ai,Bi. 

^i.n(n — k — l)(n — k — 2)...(n — 2k + 1) „ ni. i , n—1 , ^ 

A „-i2,-i) = (-1)^^ '-^ ' ^ L2n-2k-i^ k = 0,2,...,^—, 36 

2 fc! 2 

= (-l)^ ^""^"^^^""t"'^-^""'^^ 2"-^^-^ , = 0, 1, 2, ^ . (37) 
2 k! 2 



It can be verified that for n = 1, 3, the Ai, Bi which follow from Eqs. (j33p to (j35p agree with the values 
obtained by us above as well as in pp. Further, for low values of (odd) n, it can be numerically checked 
that the Ai,Bi which follow from Eqs. (j36|) and (|37|) indeed satisfy the various relations which follow by 
demanding the validity of the constraint relation ^ . 

In the next section, we will prove that our conjectures hold for all n. But first, using Eqs. (|33|) . (|36p . 
and (|37p . we predict that the following identities hold for all odd n: 

, ^•.n(n — k — l)(n — k — 2)...(n — 2k + 1) „ oi- ^ , ■. 
— :^2"-2fc-i = 1 . (38) 

fc=0 



E^i.(n — k — l)(n — k — 2)...(n — 2k) „ oi. i , , 

(-l)'^-^ ^-^ ^—^ ^2"-2'=-i = n . (39) 

k=0 

We will see that these identities naturally lead to the proof of our conjectures. 
Case II: Lame Polynomials of even order 

In this case there are n/2 + 1 number of terms in /„ and n/2 number of terms in Qn- In this case, one 
of the solutions can be written in the form (note n is an even integer) 

(n/2+l) 

fn= Yl ^fc(dn[/3(n + C2),m])2('=-i), 
k=l 

(n/2) 

gn = Vmsn[/3{n + C2),m]dn[P{n + C2),m] ^ 5fc(dn[/3(n + C2), m])^^^-^) . (40) 

k=l 

For this solution, fn,gn and hence Un,Vn satisfy the boundary condition ([IT]) . 



Another solution is given by 

n/2+1 

fn= Yl ^fc(cn[/3(?i + C2),?n])2{^-i), 
k=l 

n/2 

5„ = sn[/3(n + C2),m]cn[/3(n + C2),m]^Sfc(cn[/3(n + C2),m])2(*^-i) . (41) 

fc=i 

For this solution, fn,gn and hence Un,Vn satisfy the boundary condition ()13p . 



In the limit m = 1, both the solutions (j40p and (j4ip go over to the hyperbolic solution 

n/2+1 

fn= Y ^fc(sech[/3(7i + C2)])2('=-^), 

k=l 

n/2 

gn = tanh[/3(n + C2)]sech[/3(n + C2)] ^ 5fc(sech[/3(n + 02)]?^''"^'^ . (42) 



k=l 
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Further, in the limit m = 0, the solution ()4ip goes over to the trigonometric solution 

n/2+l 

fn= Yl Ak{cos[P{n + C2)]f('^-'K 

k=l 

n/2 

Qn = sin[/?(n + C2)] cos[/3(n + cs)] 5fc(cos[/3(n + 02)]?^^^^'^ . (43) 

Note that these four solutions are valid in case fJ-i,fJ,2 < 0. However, if fii and fi2 have opposite signs, 
say /ii < 0, //2 > 0, then the solution is given by 

n/2+1 

f = V 

^ (dn[/3(7i + C2),m])2(fc"i) ' 

ff„ = V^sn[/3(n + C2), m] ^ ^ IWra " ^^^^ 

^ (dn[/5(n + C2),mJ)^^'^ 



For this solution, fn,9n and hence Un,Vn satisfy the boundary condition (|TT|) . In the limit m = 1, this 
solution goes over to the hyperbolic solution 

n/2+1 

/„= ^ ^fc(cosh[/3(n + C2)])2('=-i), 
fc=i 

n/2+1 

5„ = tanh[/3(n + C2)] J] i?fc(cosh[/3(n + C2)])2(^-i) . (45) 

fc=2 

On substituting any of the expressions for fn,gn as given by Eqs. (HOj) to (HSj) in Eq. ([6]) we again 
obtain n + 1 equations which determine the n + 1 parameters A/j , . Some of these relations are 

^n/2+1 — -°n/2 ' 

2^n/2+l^n/2 + ^n/2 = 2Sn/2-B„/2-l , 

^n/2 + 2^n/2+l^n/2-l + 2i3„/2-B„/2-l = -B^/2-1 + 2-Bn/2^n/2-2 ) 
^n/2-1 + 2^n/2+l^n/2-2 + 2-B„/2^n/2-2 
+2A„/2^ri/2-l = 2S„/2-l-B„/2-2 + 2i?„/2-B„/2-3 ) 

Al = l, Bf + 2AiA2 = 0, 2A1A3 + Al + 2B1B2 = Bf , 

Bl + 2A2A3 + 2^1 + 252^2 = 251^2 . (46) 
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On comparing these results with the exact expressions for n = 2, 4, we conjecture the following general 
results for arbitrary even n: 

n(n^-3) („_5) (re-3)(?i-4) („_5^ 

^n/2-l — ^ > -Dn/2-2 = 2^ , 

n(n-4)(ra-5) (n - 4)(ra - 5)(n - 6) (^.^^ 

^n/2-2 = ^ ) -Dn/2-3 = ^ ) 

n/2+1 n/2 

5^ A, = l, 5]A = n. (47) 

i=\ i=l 

Further, depending on if n = 4A: + 2 or 4/c + 4 we have the following results (A; = 0, 1, 2, ...): 

Ai = -1, 5i = n,^2 = Y, i/ n = 4A: + 2, (48) 

^i = +l, Si = -n, A2 = -Y, n = 4A: + 4. (49) 
In fact looking at these general results, we conjecture the following expressions for various Ai,Bi. 

^„/2-..i = (_i)^ -(--fe-l)(--^^-^)-(--^^ + l) 2»-^^-i , . = 0, 1, 2, n/2 , (50) 



„ . ,i.(n — k — l)(n — k — 2)...(n — 2k) „ oi- i , , , . 

Bn/2^k = i-'^)'- ^1 — ^2"-2'=-i, fc = 0,l,2,...,n/2-l. (51) 

It can be verified that for n = 2,4, the Ai,Bi which followfrom Eqs. (|46p to (j49p agree with the 
values obtained independently by us above as well as in(T|. Further, for low values of (even) n, it can 
be numerically checked that the Ai,Bi which follow from Eqs. (ISOh and (15ip indeed satisfy the various 
relations, which follow by demanding the validity of the constraint relation ([6]). 

In the next section, we will prove that our conjectures hold for all n. But first, using Eqs. (|49p . (jSOp . 
and (j5ip . we predict that the following identities hold for all even n: 

n/2 

^-^k ^jn -k- l)(n -k- 2)...(n - 2k + l) ^„_2fc-i ^ ^ ^^2) 

, ^i.(n — k — — k — 2)...(n — 2k) „ n, , , , 

(-1)^^^ ^-^ ^2"-2'="i = n . (53) 

fc=0 

These identities will naturally lead to the proof of our conjectures. 

As mentioned earlier, since the field Eqs. ^ and ([5|) are invariant under /„ — )• zizfn,gn ifi'm hence 
one can trivially write down three other solutions in both odd and even n cases. 
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4 Connection to Chebyshev Polynomials 

In the previous section, we have defined two famihes fn,gn and conjectured that they lead to the solutions 
of the constraint relation Here, /„ is a polynomial, and gn equals a polynomial times an extra factor, 
and we took their common argument to be a Jacobi elliptic function. We now prove that these proposed 
functions do indeed give rise to solutions for all n. To do so, we take a closer look at the identities stated 
in Eqs. gS^I, i^, §7^1, and If we look at Eqs. dMI) and then we notice that in both the cases, 

the summation goes from k = to k = m, where n = 2m + 1 in Eq. (j38p and n = 2m in Eq. (j52|) . In 
other words, in both cases the summation goes from fc = to [n/2j. Thus Eqs. ([55j) and ([5^ change to 
the odd and even cases of 

It turns out that replacing the 2 in the numerator by other bases yields other similar identities as well, 
that can be numerically verified for small n. To simplify the notation, define 

where x € M. Our identities ()38p and (|52p say that /n(l) = 1- But before we prove this, we remark that 
we can also compute fn{x) for other values of x and low values of n. We omit writing the details down, 
but doing so leads to the following (conjectured) identities: 

/„(!) = 1 = cos(2n7r) , /„(l/2) = cos(n7r/3) , /„(0) = cos(n^/2) , (56) 

Similarly, it is possible to propose closed-form expressions for /„(— 1/2) and /„(— 1) as well. The connection 
between the terms on both sides in Eq. (j56p is made by noting that 



l = cos(27r), 1/2 = cos(7r/3), = cos(7r/2). (57) 
This leads to the following result for all x. 
Theorem 1. For any 6 G [0, tt], 

/„(cos(6l)) = cos(n6l) . (58) 
More generally, for any 2; G M, fn{x) = Tn{x), where is the nth Chebyshev polynomial of the first kind. 
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Proof: We make use of an "explicit formula" on Chebyshev polynomials: 

Simplifying the right-hand side, we obtain the desired result. Note that for x S [—1, 1], x = cos{9) for a 
unique 6 G [0,7r], and then the definition of T„ proves that /„(cos(6')) = cos{n9). □ 

It is also easy to compute that when x = 1, 1/2, or 0, Tn{x) = 1, cos(n7r/3), or cos(?i7r/2) respectively, as 
claimed in Eq. ()56p . 

Similarly, if we look at Eqs. (j39|) and ([53|) . then we observe that both these cases change to the odd 
and even cases of 

1 ^"/'J /rv / 1\ 

I )(-iy2"-^'=n. (60) 

It remains to compute the above series. More generally, we prove: 
Theorem 2. Let 

Then for all x G M, 9n{x) = xUn-i{x), where Un{x) are the Chebyshev polynomials of the second kind. 



1 L"/2J / , IN 

9n{x) = \Y.[ I ')(-i)'(2-r-^'- (61) 



Proof: We start from Eq. (j59p . On differentiating both sides, we have 

Kix) = E(-l)' 1r^;n? ("-^0(2xr-^'-^ (62) 
= -gn{x). (63) 



Simplifying, we get: 



gn{x) = - - T^^ix) = xUn-l{x), (64) 

n 



where the second equality is standard. □ 

It remains to prove the original identity. But the summation is simply gn{^), and it is well known that 
Unil) = n + 1 for all n. Hence 

Qnil) = Un-iil) = n . (65) 
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In fact it is now clear that for arbitrary n, fn and Qn (as given by Eqs. (j55p and (j6ip respectively, for 
any even or odd integer n) are simply 

fn = Tn{y), gn = {I - vY^^Un^iiy) , (66) 

where y is one of the following: 

dn[/3(n + C2),m] , cn[/3(n + C2), m] , sech[/5(n + C2)] , cos[/3(n + C2)] , 

, rof ^ cosh[/3(n + C2)]. (67) 

dn[/5(n + C2),mJ 

Note that 

fniy) + ffn(y) = T^iy) + (i - 2/')f^n-i(y) = i • (68) 

In other words, the solutions of the coupled equations are simply Chebyshev polynomials with argument 
in terms of Jacobi elliptic functions. Actually, once one realizes this, then the structure of other solutions 
is also simplified. In particular, it is known that 

Tn{cos{6)) = cos{n6) , r„(cosh(x)) = cosh(7ix) . (69) 

Thus, fn{x) = cos{n6) or cosh(nx), which makes gn{x) = sin(n0) or sinh(nx) respectively. Thus, one now 
has a better understanding of the solutions ([3T]) , ([15]) , (p9]) and ([13|) . 

5 Other Coupled Models 

We now consider the coupled Ablowitz-Ladik (AL), coupled (j)^ and coupled (j)'^ models as discussed in our 
previous publication [1] and show that all these models also admit solutions in terms of Lame polynomials 
of arbitrary order. 

5.1 Solutions of a Coupled AL Model 

As shown in our previous publication [T], in the special case when vi = 2/ii and U2 = 2/12, the coupled 
Salerno model as given by Eqs. ([1]) and ([2]) reduces to the coupled AL model with the field equations 

idUn/dt + [Un+l + Un-1 - 2m„] + (/ii | li„ ^ +^l2 \ Vn \'^)[Un+l + Un-l] = , (70) 
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idvn/dt + 



2/i2 

Vn+1 + Vn-1 Vr. 



+ (W I Un P +/U2 I Vn |^)bn+l + ^n-l] = 0. (71) 



It is then clear that aU the solutions of coupled Salerno model in terms of Lame polynomials of arbitrary 
order, are automatically the solutions of the coupled AL model and further in this case, wi = 2,uj2 = ^j^- 



5.2 Solutions of a Coupled Discrete 0^ Model 

The field equations of the coupled discrete 4>^ model discussed in our recent paper [T] are 



^((/>„+i + (fin-i - 20„) = ai(/>„ - + dil^l4>n + [cicj)^ + ecplipl + /V'^][(/>„+i + cpn-i] , (72) 
1 e 

■J^{lpn+l + - 2V'„) = a2tpn - ^2^^ + dfyn + [C2lpn + l^'Pi + 2/</'n'0n] [V'n+l + V'n-l] ■ (73) 

As shown in [1], solutions to these coupled equations in terms of Lame polynomials of order one and 
two are obtained in case 



<tl + ^l = \l^, (74) 



and further if 



ci = C2 = / = I , h = b2 = -d, ai = a2, cih^b^ = 1 , ai + ^ = ^1— . (75) 
2 ft^ \/h^ci 

It is then clear that the solutions (|26p to (|'28p and (|40|) to (|42|) in terms of Lame polynomials of arbitrary 
order obtained in the case of the coupled Salerno model will also be the solutions of this coupled model. 

5.3 Solutions for a Coupled Discrete cf)'^ Model 

In our recent publication |T] we considered the following coupled (j)^ model 

^{(t>n+l + 4>n-l - 20„) - 2ai^n " [2pl4>l + li^l][<Pn+l + 4>n-l] = , (76) 
^(^„+l + - 2^n) - 2a2V'n " {2^2i^l + 7<^n][V'n+l + i'n-l] = , (77) 

and showed that solutions to these coupled equations can be obtained in terms of Lame polynomials of 
order one and two provided 
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and further if 

2/3i = 2/32 = 7, ai = a2 = -^. (79) 

It is then clear that the solutions (|26p to (|29p and (|4Up to (|43p in terms of Lame polynomials (or corre- 
sponding hyperbolic and trigonometric polynomials) of arbitrary order obtained in the case of the coupled 
Salerno model will also be the solutions of this coupled model. 

6 Summary 

In this paper we have shown that for a number of coupled discrete models, e.g., coupled Salerno, coupled 
Ablowitz-Ladik, coupled 0^, coupled (j)^, there are solutions in terms of Lame polynomials of arbitrary 
order while the uncoupled equations do not admit solutions in terms of Lame polynomials of order two and 
higher. In particular, we showed that the Lame polynomials can be reexpressed as Chebyshev polynomials 
of the relevant Jacobi elliptic function. Many of these solutions are relevant to physical contexts ranging 
from ferroelectric materials [HJ [9l [10] to multiferroics O [U [7] to a variety of models in field theory [TTl [15] 
in addition to various discrete contexts [HI [131 [H] • 

It is important to emphasize that the stability of the various solutions found here remains an open issue 
to be explored numerically, particularly since the soliton solutions obtained above are of arbitrary width. 
Besides, a comparative study of the scattering of solitons of different discrete models is an important issue 
with these static solutions boosted with a certain velocity. Similarly, the Peierls-Nabarro (discreteness) 
barrier for the solutions remains to be explored. However, since all our solutions have discrete translational 
invariance (i.e., they are valid for arbitrary C2), it is likely that for all our solutions the Peierls-Nabarro 
barrier may be zero. However, this issue needs to be explored carefully and we intend to do so in the near 
future. 
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